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Context. The interpretation of the X-ray spectra of X-ray binaries during their hard states requires a hot, optically 
thin medium. There are several accretion disc models in the literature that account for this aspect. However, none is 
designed to simultaneously explain the presence of powerful jets detected during these states. 

Aims. A new quasi-keplerian hot accretion disc solution, a Jet Emitting disc (JED hereafter), which is part of a global 
disc-jet MHD structure producing stationary super-alfvenic ejection, is investigated here. Its radiative and energetic 
properties are then compared to the observational constraints found in Cygnus X-1. 

Methods. We solve the disc energy equation by balancing the local heating term with advection and cooling by syn- 
chrotron, bremsstrahlung and Comptonization processes. The heating term, disc density, accretion velocity and magnetic 
field amplitude were taken from published self-similar models of accretion-ejection structures. Both optically thin and 
thick regimes are considered in a one temperature gas supported disc. 

Results. Three branches of solutions are found possible at a given radius but we investigate only the hot, optically thin 
and geometrically slim solutions. These solutions give simultaneously, and in a consistent way, the radiative and ener- 
getics properties of the disc-jet system. They are able to very well reproduce the global accretion-ejection properties of 
Cygnus X-1, namely its X-ray spectral emission, jet power and jet velocity. About half of the released accretion power 
is used to produce two mildly relativistic {v/c ~ 0.5) jets and for a luminosity of the order of 1% of the Eddington 
luminosity, the JED temprature and optical depth are close to that observed in the hard state Cygnus X-1. 
Conclusions. The accretion and ejection properties of JEDs are in agreement with the observations of the prototypical 
black hole binary Cygnus X-1. The JED solutions are likely to be relevant to the whole class of microquasars. 
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1. Introduction 

Powerful ejections are commonly observed in galactic black 
holes. X-ray binaries show perman ent jet radio emission 



1977). In this model, accretion serves only to bring in and 
maintain the magnetic field in the black hole ergosphere. 

The fact that the jet power observed in some X- 



during the so-called hard states ([Stirling et al 



2001 
Remillard fc McClintock[ |2006[ [Done et al.[ |2007| ) and"iTl 



ray binary syste ms (the so-called microquasars Mirabel & 
Rodriguez||1998 1 seems remarkably independent of the na- 



perluminai sporadic e jcctiWhf intcrmedia teilmng states tiF^^HEi compact objects namely white dwarf neutron 



( Mirabel fc Rodrigu ez 
Fomalont et al.||2001 



1994; Hannikainen 
Fei^deFet al.||2004[ P 



et al. 2001 



Several models have been put forward to explain the 
formation of these astrophysical jets. They can be roughly 
classified in two categories both relying on the existence of 
an organized large scale magnetic field. The first one taps 
th e accretion energ y res ervoir of the disc, as first proposed 
by Lovelace (|1976^ and Blandford ( 1976 ). Since the semi- 
nal paper of Blandford & Payne ( 1982p , it is known that 
jets can indeed be both accelerated ana coUimated by the 
magnetic field threading the disc. In the second category 
of models, jets are supposed to mainly tap the rotational 
energy of the black hole thanks to the so-called Blandford- 
Znajek process (hereafter BZ process, Blandford fc Znajek] 



^ A tentative summary of our knowledge of 
accretion-ejection ph enomen a can be found here 



http://www.issibern.ch/teams/proaccretion/ 



star or black hole (e.g. Kording et al. 2006 Kording et al. 



2008 Tudose et al. 2009) argues for an ejection process 



independent of the central object. This seems to be con- 
firmed by the lack of correlation between jet powe r or jet 



velocities and measured values of black hole spins ( Fender 



et al. 2010 ). Although all these evidences do not constitute 
a definitive proof, it severely weakens the scenario where 
powerful jets would be powered by the BZ process. 

There have been numerous 2D and 3D relativistic niag- 
netohydrodynamic (hereafter MHD) numeric al simulations 
of accretion fiows around black holes (e.g. |McKinney 
Gammie 2004; Komissarov ct al. 2007; Hawley 



& 
Krolilc 



2006;iBeckwith et al., 2008a. .McKinney fc Blandford,2009p 
Starting from a geornetrically thick donut pervaded by its 
own weak poloidal or toroidal magnetic field, simulations 
end up with accretion onto the black hole with outflow- 
ing material along opened field lines. Accretion is due to 
the onset of the magneto-rotational instability or MRI (see 
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Balbus 2003 and references therein) that leads to the for- 
mation oFmi almost geometrically thin disc in the inner- 
most radii. A Poynting-dominated funnel jet is formed in 
the black hole polar regions in good agreement with the 
BZ process. This electromagnetic (no mass) funnel wind is 
surrounded by a mass-dominated outflow concentrated in a 
hollow and relatively thin cone at roughly 40 degrees incli- 
nation. Ejection is both thermally and magnetically driven 
and is only prese nt if the initial magneti c field had a ver- 



tical component (Beckwith et al. 2008a). Although these 



simulations look very promising, it is still not clear how 
they compare to the observed large scale jets. Indeed, it is 
not clear yet whether such flows will be self-confined fur- 
ther out and if the mass and energy fiuxes are in agreement 
with observational constraints. 



On the other hand, following the Blandford fc Payne| 
( 1982 1 pioneering work, it has been realized that jets 
could extract a significant fraction of the underlying disc 
angular momentum and accretion power. Under such 
circumstances, jets would not be a mere epiphenomenon 
of accretion but possibly its main driver. The whole disc 
structure should then be revisited by taking into account 
the disc-jets interrelation. This was done within a self- 
similar Ansatz, allowing to solve the full set of dynamical 
MHD equations without neglecting any term (Ferreira &' 
Pelletier[ [19951 |Ferreira[ p^97l [Casse fc Ferreirat 2000a|bi 



Ferreira fc Casse |2004[ ). It was found that the structure of 
these Jet Emitting Discs (hereafter JED) is indeed deeply 
modified, with dynamical properties significantly different 
from those of a standard accretion disc or an ADAF. As 
a consequence, the radiative properties of a JED may be 
also very different. Noticeably, the main results of these 
steady-state self-similar studies have been confirmed by 
2D numerical MHD simu l ations done with the VA C code 
dCasse fc Keppens[ |2002[ I2004[ ) , t he FLASH code dZanni 



et al.[|2007p and the PniTOrode ( [Tzeferacos et al.[|2009p ." 



The purpose of this paper is to detail the dynamical, 
energetic and radiative properties of Jet Emitting Discs 
around compact objects and to compare them with those 
estimated from the observations of the prototype of black 
hole binaries: Cygnus X-1. 

The structure of the paper is as follows. We start by 
recalling the main characteristics of a JED in Sect. |2] We 
also present in Sect.|3]the method used to compute the JED 
thermal equilibrium. We then compare, in Sect.|4]the model 
expectations with the observational constraints known in 
Cyg X-1. Concluding remarks are discussed in Sect.[5J 



2. Jet Emitting Discs 



Jet Emitting Disc s have been originaly studied by [Ferreira] 
& Pelletier ( 1993 1 in their Magnetized Accretion Ejection 
Structures model (MAES herafter). This model has been 
developed so as to treat consistently both the accretion 
disk and the jet it generates. The idea is the same as in 
earlier studies of magneto-c entrifugally launched disk winds 
( Blandford fc Payne[ [T982| . However, in the MAES, the so- 
lution starts from the midplane of the resistive MHD disk 
and evolves outwards in the ideal MHD wind/jet. This dif- 
fers drastically from other studies where the disk was only 
treated as a boundary condition, hence forbidding any pre- 
cise quantication of the effect of the MHD wind on the disk. 



It would be lengthy but also irrelevant to present the 
MAES model in great detail in this paper. Many papers 
have dealt with the subject, from both analytical and nu- 
merical point of views, and we refer the reader to these pa- 
pers for further details (e.g. Ferreira|1997 Casse & Ferreira 
2000a|b[ JFerreira[[2002l [Ferreira fc Casse||2 004; Casse fc 
Keppensl'2002[ [2004[ [Zanni et al.[[2007|[Combet fc Ferreira 



^OOSj). Instead, we give hereafter the few key assumptions 
and elements of the model that are important to our work. 



2.1. Main assumptions 

The goal of these studies was to find out the conditions 
for steady state jet formation and to relate jet properties 
(mass loss, power, kinematics, coUimation) to the under- 
lying disc properties. The MAES calculations were done 
under several assumptions that are discussed below: (1) a 
self-similar Ansatz; (2) an alpha-prescription for the trans- 
port coefficients; (3) the presence of a large scale, organized 
vertical B^ field. 

(1) Self-similarity enables to solve the full set of MHD 
equations without ignoring any term. The only simplifica- 
tion that must be made is in the energy equation, since mi- 
crophysics (within cooling terms) is not self-similar. Thus, 
given some approximation (i.e. isothermal or adiabatic 
magnetic surfaces), a computed solution is a real solution 
of the full set of equations. In our case here, such a solution 
is obtained by crossing the MHD critical points of the fiow. 
On the other hand, self-similarity introduces una voidable 
biase s in the jet coUimation properties (see e.g. Ferreira 



1997 for more details) 

(2) The disc is assumed to be turbulent so that, within 
our mean field approach, some prescription must be made 



to mimic turbulence. Following Shakura & Sunyaev (19731), 
an alpha prescription for the transport coefficients, namely 
viscosity and magnetic diffusivity, is employed. It has been 
shown recently that MRI gives rise to a turbulent magnetic 
transport that seems to behave like a resistivity, with an 
effective ma gnetic Prandtl number of ord er unity ( Guan fc 
Gammie, 2009 Lesur & Longaretti 2009). It turns out that 
JED solutions are found to be close t o the marginal stabil 
ity limit, but always MRI unstable (Ferreira & Pelletier 



1995 Ferreira 19971. Note that there is a whole range of 



other possible MHD instabilities when the fields are close to 



equipartition (Keppens et al. 



2007). Whether magnetically driven turbulence can indeec. 



2002 Blokland et al. 2005 



be descri bed with local alpha p rescriptions is still an open 
question (Beckwith et al. 2009). Full 3D global simulations 



are thus required to assess this assumption, but the huge 
span in space-time scales are such that this is still out of 
reach of modern computers. 

(3) The last assumption, namely the presence of a large 
scale magnetic field in the inner regions of accr etion flows, is 
still an open and highly debated question (e.g. [Lubow et al. 



1994a|bl [Cao fc Spruit||2002l [Bisnovatyi-Kogan fc Lovelace 



|2007[ [Rothstein fc Lov elace"2008[) . 

It has been recently argued However that fields could be 
advected inward in a standard accretion disc thanks to the 



quenching of the MRI at the disc upper layers ( Rothstein 



& Lovelace 2008 Lovelace et al. 2009). This picture has 



been given some support with recent 2D MHD simula- 



tions of accretion-ejection structures (Murphy et al. 2010 
GRMHD simulations also suggest that a vertical large sea 
magnetic field is indeed needed in order to produce power- 
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ful jets ( Beckwith et al. 2008b ) . Thus while the controversy 
is certainly not closed, tne assumption of a large scale mag- 
netic field threading the disc appears quite plausible. 



where we introduced two dimensionless parameters, the 
sonic Mach number rrig — Uq/cs and the disc magnetiza- 
tiofji = fio'B^jP. 



2.2. Physical properties 

Throughout the paper, we use R for the cylindrical radius, 
M for the black hole mass and Ma{R) for the disc accretion 
rate at the radius R. The dimensionless radius will be r = 
R/Rg where Rg — GMjf?, mass to = M/Mq and disc 
accretion rate m 
luminosity. 



MaC^/LEdd with Lsdd the Eddington 



2.2.1. Radial profiles 

Since a JED undergoes mass loss, the disc accretion rate is 
written as 



Ma{R) = Ma,. 



R 

^out 



(1) 



where ^ measures the local disc ejection efficiency (Ferreira 



& Pelletier 1993) and Mout the accretion rate at the disc 
outer radius Rout- A standard accretion disc is recovered 
when 'C = 0- The ejection efficienc y ^ is equivalent to the p 



expon ent used in ADIOS models (Blandford & Begelman 
19991. But, in strong contrast to the latter, it is not as 



sumed but computed as function of the d isc parameters as 
a trans- Alfvenic regularity co ndition ( see Ferreira||l997 for 



more details). For isothermal ('Ferreira 1997) and adiabatic 
(Casse & Ferreira 2000aj magnetic surfaces, ^ is found to 
vary in a narrow interval around 0.01. When the disc upper 
layers are heated, mass loss can be enhanced with ^ reach - 
ing typically 0.1 or slightly more (Casse & Ferreira 2000b). 



Discs with larger mass losses do not give rise to stationary 
super- Alfvenic jets. 

The disc vertical equilibrium depends on the balance be- 
tween, on one side, gravity and magnetic compression due 
to both radial and toroidal field components and, on the 
other side, the vertical pressure (gas -I- radiation) gradient. 
We nevertheless use as a proxy for the real scale height the 
hydrostatic value, namely P ~ po^\H^ evaluated at the 
disc midplane, where P = Pgas + Prad'- this greatly sim- 
plifies our expressio ns and produces an overestimation of 
at most a factor 2 (Ferreira & Pelletier 1995). With this 
caveat in mind, the disc aspect ratio in a gas supported 
disc is e = H/R — c^/vk, namely the ratio of the isother- 
mal sound speed to the Keplerian speed vk = ^kR- Given 
these definitions, the radial profiles of the accretion veloc- 
ity Uo, particle density n — po/rrip, gas pressure Pgas and 
vertical magnetic field B^ in the JED midplane are given 

by 



-Ur = rrisCs — rtis e Vl^R 

Ma{R) 



P — 

^ gas — 



AnrnpflK R^TTT-s^"^ 

2 MaiR)nK 

PoCs = 



B, = 




p \ ■ J PoMa{R)VLK 

TO, / V 47r_R 



(2) 
(3) 

(4) 
(5) 



2.2.2. JED parameter space 

The previous expressions ([2])- ([5]) are common to all power 
law, nearly Keplerian gas supported disc solutions. What 
distinguishes JEDs from other models are the values taken 
by the parameters ^, to^, p and e. 

— The value of e is fixed by the resolution of the JED 
energy equation (see Sect.js]). 

— p and ^ are provided by the resolution of the full set of 
MHD equations when properly dealing with the regu- 
larity conditions at the slow magnetosonic and Alfvenic 
critical points, respectively. We already discussed the 
values achieved for ^. For isothermal or adiabatic mag- 
netic surfaces, steady-state ejection is possibl e only for a 
field smaller than but close to equipartition ( Ferreira fc 



Pelletier! [T9951 [Ferreira] |1997[ [Casse fc Ferreira[ |2000ap 
More precisely, the range of allowed magnetizations m 
isothermal models goes from 0.3 to 0.8, with an aver- 
age value aro und 0.6 (see examples of JED parameter 
values in Fig[A.l[). This can be understood in the fol- 



lowing way. The magnetic field is vertically pinching the 
accretion disc so that a quasi-static vertical equilibrium 
is obtained only thanks to the kinetic pressure support. 
As a consequence, the field cannot be too strong. But 
on the other hand, the field must be strong enough to 
accelerate efficiently the plasma at the disc surface so 
that the slow-magnetosonic point is crossed smoothly. 
In a JED, the sonic Mach number can be rewritten un- 
der the form 



2qp 



(6) 



where a^e denotes the effect of standard transport and 
2qp the specific contribution of the magnetic torque. 
The magnetization p = B^ / p^P measures the strength 
of the magnetic field in the disk and q = poJrh/ B^ 
is the normalized radial electric current density flowing 
at the disk midplane. This last parameter measures the 
magnetic shear as it provides an estimate of the toroidal 
magnetic field component at the disk surface, namely 
B+ c -qB^. 

In a SAD, tti^ = aye, with typically q;„ = 10^" and 
e — H/R < 1, so that the accretion velocity is largely 
subsonic. However, the situation in a JED is very differ- 
ent since the presence of powerful jets strongly modified 
the disc angular momentum. A common feature of all so- 
lutions found is that the jet torque always domi nates the 
local turbu l ent to rque, as initially proposed by [Pelletier[ 
fc Pudritz[ ( [1992[ ). This is because steady-state MAES 
solutions are generally found close to equipartition (i.e. 
p ~ 1 as explained before) and g also of the or der of 
unity (Ferreira fc Pelletier 1995 Ferreira 1997). As a 



consequence JED solutions present high accretion veloc- 
ities with TOs ~ 1. And this i s ver ified by all isothermal 
and adiabatic solutions (Fig. A.l). 



In order to simplify the analysis carried out in this pa- 
per, we will fix hereafter rris = 1 and p =^ 0.5. Putting 



^ In a gas supported disc, p = 2//3 where /3 is the well known 
plasma beta parameter. 
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numbers, one gets: 

Pgas — 1-5 lO^*" mm~^r~^/'^ erg.cm"^ 
B^ ~ 4.4 10® ml/2 TO-i/2 ^-5/4 Q 



(7) 
(8) 
(9) 



Blandford & Payne (19821 scaling for the 



Note that the 

magnetic field is recovered as we used ^ -C 1. 



2.2.3. Global energy budget 

The gfobal energy budget in an accretion-ejection structure 
is 



Pa 



PjED 



p. 



jets 



Pn 



dv 



p. 



jets 



(10) 



where Pjets is the total power feeding the jets and Pjed 
is the power dissipated within the JED. More precisely, 
Pjed is the sum of the radiated Prad and advected Padv 
power. The disc luminosity is therefore the released accre- 
tion power minus what is advected onto the black hole and 
what feeds the jets. 

Given the expression (fTl) of the accretion rate, the total 
accretion power released in a quasi-Keplerian accretion disc 
writeC] 



p — 

^ ace 



GMMa{R,n) GMMajRo 

GMMa,out 



where 



(1-5) 



GMMg^OUt 



Rii 



Ri- 



Ro 



^out 

{1-9) 



2i?o«t 



^in 



Rn 



Ro 



(11) 



(12) 



namely Paee = V^ ^Edd with an accretion efficiency 
r/ = (1 — g)/2ri„. The usual expression for the non rela- 
tivistic accretion efficiency is recovered for ^ = 0. In the 
following, we will use ri„ = 6 for a Schwarzschild black hole. 

The power Padv that is advected onto the central object 
along with the accreting material scales as Padv oc e'^Paec- It 
is thus negligible in a thin disc. Pjets is the sum of the total 
(kinetic, potential and internal) energy flux advected by the 
outflowing plasma and the MHD Poynting flux leaving the 
disc surfaces. Since we are mainly interested in powerful 
jets, the latter contribution dominates i.e. 



Pjets — Pmhd 



n,R 



^J■o 



'-2t: RdR 



(13) 



where Bf is the toroidal component of the magnetic field at 
the disc surface and fi* is the rotation speed of the magnetic 



This expression neglects any input of energy from the cen- 
tral object at the disc inner e dge, through turbulence or MHD 
Poynting flux (see for instance Casse & Ferreira|2000b I . 



surface^ Using the disc angular momentum conservation, 
one gets (see Appendix [A| 



P. 



jets 



bPae 



(14) 



where 6, the jet power parameter, is known for a given MHD 
solution. Note that in the above derivation, it has been as- 
sumed that the fraction of the released energy that is trans- 
ferred to the jets is a constant throughout the disc. That 
is not necessarily true. Indeed, self-similar models (where b 
is by definition a true constant) published in the literature 
provide b ranging from ^0.5 to almost unity. It turns out 
that this value depends on the disc aspect ratio but its ac- 
tual dependence remains to be looked for. As a first step 
and for the sake of clarity in this paper, we will use it as 
an independent parameter. 

Consequently, given the above expressions, the total 
JED power Pjed and JED luminosity Prad are given by 



P 



JED 



^acc ^jets 



P 



JED 



P. 



dv 



(1 - b)Pacc 
= (1 - b)Pa. 



idv 



(15) 
(16) 



The jet power parameter b is thus a crucial parameter as 
it controls the power sharing between the jet and the disc. 
Sinc e it is generally found in the range 0.5-0.99 (see Fig. 
A.l I, JEDs considered here are engines converting accretion 



power into ejection power with a high efficiency. To get the 
fraction of the power Prad that is actually radiated away, we 
need however to compute the JED thermal balance. This 
is explained in the next section. 

3. Disc thermal balance 

The equation governing the internal energy of the accretion 
fiow writes 



Qheat = qrad + PV.U + V[/.U 



(17) 



qadv 



where qheat is the heating power density, qrad the sum of all 
radiative cooling terms, P the total plasma pressure, u the 
fiow velocity, U its internal energy and qadv the advection 
term. These different ter ms a re explicited below. Equation 
p^l is equivalent to Eq. ( 16 ) but expressed locally in the 
JED. 



Note that in Eq. (17) we neglect any turbulent en- 



ergy transport (such as convection for instance) and as- 
sume a Maxwellian one temperature plasma. This last 
point could be a crude assumption at low accretion rates 
(to < lO^"' — 10— 4). However, since we are mainly interested 
in this paper to apply our calculations to Cyg X-1, whose 
accretion rate is of the order of a few % of the Eddington 
rate , this assumption is expec ted to be quite reasonable (see 
e.g. [Malzac fc Belmont||2009[ ). 



3.1. Heating and advection terms 



Equation ( 15 ) provides us directly the radial JED volumet- 



ric heating rate in a disc ring of height 2H, radius R and 



'^ A steady-state jet is a bunch of axisymetric magnetic sur- 
faces nested around each other and in differential rotation. Each 
magnetic surface is put into rotation by the underlying disc, so 
that SI* ~ r2^(r), where il -I- (r) is the angular rotation at the 
disc surface. Note that this velocity can be significantly smaller 
than that at the disc midplane. 
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extent dR: 
qheat{R) = 



1 dPjED 

AttRH dR 

(1 - b) dPac 



AnRH dR 



= (1-6) 



GMMa 

8ttHR^ 



(18) 
(19) 



-1.5x10^^(1-6) — ^-r erg.s-i.cm-3. (20) 

On the other hand, the advection term is identified as 

1 



Qadv = PV • U + VC/.U 



7-1 



[7PV • u + u .VP] (21) 



where 7 = 5/3 is the gas adiabatic index. Assuming that 
the flow poloidal velocity is dominated by its radial compo- 
nent Ur (provided by Eq. pi) and after some algebra, this 
expression becomes 



m„c^m^e^n ( 9 hie 
qadv = ^tt^tI ( 4,^^^ - 1 



i?5/2 



-3 X 10^ 



SlnP 

. ^ In (_ 
o 4 , 4— 1 I erg.s \cm 



(22) 



The comparison between Eq. ( 20 ) and ( 22 ) shows that ad- 
vection becomes important only for large aspect ratios, 
namely e > 0.1. Note also that, as already examplified by 



Yuan (20011; Yuan et al. (20031 in ADAF, the advection 



term may change sign and become negative, pl ayin g then 
the role of a heating term for the flow. Equation ( 22 ) shows 



that this should occur when the exponent in the radial de- 
pendency of e becomes larger than 1/4. 

3.2. Radiative cooling term 



We strictly follow Esin et al. ( 1996 ) for the computation of 
the different radiative cooling rates. In the optically thin 
regime, the total cooling rate is simply 



Qrad QC^syjic ' QC^brem 



(23) 



where qc,sync and qc,brem correspond respectively 
to the synchrotron-|- comptonized synchrotron and 
bremsstrahlung-|-comptonized bremsstrahlung cooling 
rates. Since JED solutions may also be optically thick, 
we use the generaliz ed radiative cooling formula used by 
Narayan fc Yi (1995), and adapted from the initial work of 
Hubenyl 



(|1990l, that bridges the optically thin and thick 



cases: 



Qrad 



l.bT + VS- 



[AaT^/H){qc^sync 



QC,bre7n ) 



(24) 



^ Tabs is the disc optical depth in the ver- 
Tes being the Thomson optical depth and 

4\ 



where r = Tg^ 

tical direction 

Tabs = {H/'^cFBT'^){qc,sync + qc,brem) the Optical depth for 

absorption. At high energies, pair creation provides another 

means of cooling. But an a posteriori check showed that this 

effect remains always small and has thus been neglected. 

3.3. Existence of tliree solutions 

Figure (IT]) shows the different heating, cooling and advec- 
tion terms as function of the temperature for a given black 




6 8 

log T (K) 



10 



12 



Fig. 1. Thermal balance of a Jet Emitting Disc at r = 6 
for m =0.001 and m = 10. The underlying JED structure 
has a jet power parameter b = 0.5 and an ejection effi- 
ciency ^ = 0.1. The dotted line is the heating term Qheat, 
the dot-dot-dot-dashed line the radiative cooling term g^od 
and the long dashed line is the advective term qadv The 
comptonized Bremsstrahlung qc,brem (short dashed) and 
the comptonized synchrotron qc.sync (dot-dashed) radia- 
tive terms are also shown. The soli d lin e corresponds to 
qrad + qadv and the solutions of Eq. ^7n are indicated by 
crosses. 



hole mass, JED radius and accretion rate. In a gas sup- 
ported one-temperature disc, the midplane temperature is 
given by 



T 



IliipC- O _ 

-e r 



2k 



1 = 5.4x lO^^e^i 



K 



(25) 



The resolution of the thermal equilibrium (Eq. 17) gives 
three branches of solutions, indicated by crosses on the fig- 
ure The hottest one corresponds to the JED solution dis- 
cussed in this paper (a more detailed discussion of the other 
solutions is postponed to a future work). It is a hot, opti- 
cally thin, thermally and viscously stable disc solution that 
behaves radiati vely in a simi lar way t han the LHAF so lu- 
tions studied by Yuan ( 2001| ) (see also Yuan et aL||2006| for 
a discussion on hot one-temprature accretion flows). The 
major difference concerns the underlying dynamics of our 
JED solutions, that self-consistently include the presence 
of powerful self-collimated jets. 

Fig. [2] also indicates the different domains of JED solu- 
tions in the to -r plane and for the two extreme values of 
the jet power parameter 6, 0.5 and 1. In the latter case, a 
cold branch exists in the whole JED whatever the value 
of the accretion rate while the hot branch is only valid 
below a critical accretion rate that depends on the JED 
radius (solid line in Fig. [2]). For b = 0.5, the parameter 
space for the cold branch is now limited to accretion rates 
> lO^'' — 10^^ (bottom dashed fine) while the upper limit 
for the hot branch has shifted to slightly larger accretion 
rate (top dashed line) compared to 6 = 1, reaching to of a 
few for JED radius larger than '-^lOO Rg. A more detailed 
discussion of the JED solutions is presented in FIO. 
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Fig. 2. JED solution domain in the m -r plane. The solid 
line, corresponding to a jet power parameter 6=1, sepa- 
rates this plane in 2 regions: the top one where only the 
cold solution is valid and the bottom one where the 3 so- 
lutions exist. The region where only a hot solution is valid 
exists at very low accretion rate not shown on this figure. 
The limits for b = 0.5 are plotted in dashed lines. In this 
case the plane is divided in 3 parts with a region, at low 
accretion rate, where only the hot solution exists. 



4. Comparison to Cygnus X-1 

4.1. Observational constraints 

The energetics of the jets and the X-ray corona of C ygnus 



X-1 have been investigated recently by Malzac et al. ( [2009 
herafter M09). It is one of the best observed black-hole X 
ray binary composed of a massive star and a black hole of 
about 10 solar masses. It is also one of the rare BHB for 



which a jet was directly observed in the hard state (Stirling 



et al.[ [2001 1 . But more importantly deep radio observations 
of the held of CygX-1 resulted in the discovery of a sh ell-like 
structure which is aligned with the resolved radio jet (Gallo 



et al. 2005). Assuming that this large-scale structure was 
inflated by the inner radio jet, its kinetic power was esti- 
mated to be o f the order of the bolometric X-ray lum inosity 
of the binary (iGallo et al.l [20051 [Russell et al.| (20071) . With 



these constraints it is possible to deduce a rough estimate 
of the ratio j = Pjets/Lh of the total jet kinetic power to 
the typical X-ray luminosity in the hard state (M09) 



0.45 < j = -^ < 1.5 



(26) 



Like most X-ray binaries, the bolometric luminosity of Cyg 
X-1 does not change dramatically during the transition to 
the soft state (see e.g. [fVontera et al.||200i jZdziarski et al 



2002[ [Malzac et al.|2066[ [Wilms et ar.||2006| . T'he observed 
luminosity jumps reac h at most a factor Lg/Lh ~ 3 — 4 



( Zdziarski et al. 
likely triggered 



2002|. The transition to soft state being 

the 



Dy an increase in mass accretion rate 
ratio of the soft to hard radiative efficiencies can be safely 
upper limited by 



A 






< 4 



(27) 



It is noteworthy that there are also several estimates 
of the terminal jet velocity of Cyg X-1. Indeed, based on 
the ab sence of detection of the counter jet, [Stirhng et al. 
'2001) give a lower limit on the bulk velocity of the ra- 
10 jet of Voo/c > 0.6. Similar considerations and the lack 
of response of the rad io emission on short time-scales led 
Gleissner et al. (2004) to constrain the jet velocity in the 
range 0.4 < Voo/c < 0.7. [Ibragimov et al.[ ( [2007 1 find a 
similar result, 0.3 < Voo/c < 0.5, by modeling the super- 
orbital modulation observed in X-ray and radio bands. This 
gives another independent constraint that will be very use- 
ful while confronting models to observations. 



4.2. Model comparison 

These observational constraints on j and A can be easily 
translated into constraints on our JED parameters b and 
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Fig. 4. Average terminal jet velocity versus the disc ejec- 
tion efficiency ^ in a JED, with Ma oc R^ . The grey and 
dark grey areas correspond to the regions that agree with 
the observational constraints reported in Fig. ([s]). 



^. From Eqs. ( 14 ) and ( 16 1 and assuming that advection is 
negligible (which is vermed for JEDs at m of a few % of 
the Eddington accretion rate, as expected for Cyg X-1) we 
obtain 



Fig. 5. Contour plots of the JED temperature (solid lines) 
and optical depth (dashed lines) in the 5, — h plane at a 
distance of 10 Rq from a 10 Mq black hole. The ratio of 
the outer to inner JED radii is taken equal to 100 and the 
total accretion rate for the computation of the disc thermal 
balance is fixed to 1% of LEdd, i-e. the typical luminosity 
of Cyg X-1. The grey and dark grey areas correspond to 
the regions that agree with the observational constraints 
reported in Fig. ([s]). 



J = 



1-6 



(28) 



The soft state is free of ejection (i.e. ^ = 0) so that 
Ls — rjsms LEdd, with rj^ ~ l/2ri„. The hard state is char- 
acterized by Lh = {l-b)T]hThh LEdd, with rjh = (l-ff)/2ri„ 



where g in the hard state is defined in Eq. ( 12 ) and then 
depends on ^ and the JED radial extent RoutjRin- Within 
this framework, we obtain 



A 



1 



{l-h){l-9) 



(29) 



The left panel of Fig. (pi) displays the domain in the 
observed parameter space j — X allowed by the obser- 
vations. Conservative values are shown in grey whereas 
most plausible ones, i.e. j — 1 and A ~ 3 (jMalzac et al. 



20091, are in dark grey. The right panel shows the same 



constraints but translated into our JED parameter space 
£, — b, assuming Rout/Rin = 100. This value is com parable 
to what could be deduced from observations ( e.g. IFabian 



et al ||1989[ [Young et al.||2001| [Torres et al.||2005| ) if the 
standard accretion disc inner edge is dire ctly assimilated 
to the JED outer radius as proposed in [Ferreira et al.| 
(2006). Note that this figure depends only slightly on 
this parameter with a smooth increase of the ^ upper 
limit with decreasing Rout/Rin- We can thus safely 
conclude that ^ has to be lower than ~0.1 and b rang- 
ing between 0.4 and 0.6 to be consistent with Cyg X-1 data. 

Now, the jet terminal velocity can be also estimated. 
Since the jets of Cyg X-1 are not highly relativistic, we 
can assume that the jet power Pjets is asymptotically dom- 
inated by the kinetic energy of the ejected matter. This 
provides an average terminal bulk Lorcntz factor 



(7oo) ^ 1 + 



P. 



jets 



M.etsC'^ 



(30) 



where the total mass loss rate Mjets is such that 



/ 



M.ets 



Ma{Rout) 



Ri: 



^ou 



^In 



Ro 



^in 



(31) 



the last equality hol ding only for small ^. Combining this 
expression with Eq. (14), we finally obtain 



(7o 



1 






b (1-5) 



2r,, 



f 



(32) 



The range of average jet velocities Voo/c consistent with the 
values of b and ^ needed to agree with Cyg X-1 observations 
can then be computed from this last expression. Figure 
^ shows Voo/c as function of the disc ejection efficiency 
t~ The grey area corresponds to the conservative values 
used in Fig. (l3]), the dark grey area to adopted ones. These 
last ones provide a range in Voo/c between 0.3 and 0.6 in 
good agreement with observations. This is remarkable as it 
arises from another, independent observational constraint. 

Our treatment of the disc thermal balance allows us to 
give an estimate of the expected temperature and optical 
depth within a JED. Contours of these two quantities at 
a radius of 10 Rq are shown in Fig. (Isl), the solid lines 
representing the temperature while the dashed lines the 
optical depth. The black hole mass is taken equal to 10 
Mq and the total luminosity is about 1% of the Eddington 
luminosity. This corresponds to the typical luminosity of 
Cyg X-1 close to the hard-to-soft state transition. The ratio 
of the outer to inner JED radii is fixed to 100 but Fig. ^ 
does not strongly depend on this parameter. Interestingly, 
the temperature is of the order of 120-140 keV and the 
optical depth r ~ 0.7 — 0.9, which are not so far, given 
the simplicity of our thermal balance computation, from 
the values dedu ced from sophisticated fits of Cyg X-1 in 
the hard state ( Salvo et al. 2001 Gierlinski et al. 1997 
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Frontera et al.[ [200Ti ribraRimov et al.| [20051 |Bel et al.[ 



2006 Wilms et al. 2006). 



y 



At 



kT^ 



5. Concluding remarks 

Jet Emitting Discs (JED) can be much less radiative than 
Standard Accretion Discs fed at the same accretion rate. 
Indeed, JEDs redistribute the released accretion power into 
three outputs: radiation, advection onto the central object 
and jet power. Using published self-similar calculations of 
JEDs, which provided us with the relevant disc parame- 
ters, we were able to solve the disc thermal balance around 
compact objects in a consistent way. 



As in other studies (e.g. Chen et al. 19951, we found 
that there are usually three branches of solutions at 
a given radius: (1) a cold solution, optically thick and 
geometrically thin; (2) an intermediate solution and (3) a 
hot solution, optically thin and geometrically thicker. Only 
the cold and hot solutions are both thermally and viscously 
stable. In this paper, we focussed on the properties of the 
hot solution as it corresponds to the physical situation en- 
visioned during the hard states of X-ray binaries, which are 
known for harboring jets. We confronted our JED model 
to observations by comparing the disc X-ray emission, 
jet power and jet velocity with those derived in Cygnus X-1. 

We found that these hot JED solutions reproduce 
well both dynamical and radiative properties of Cygnus 
X-1. These solutions are characterized by a jet power 
parameter b ^ 0.5 and an ejection efficiency £, < 0.1. 
The first condition arises from the observational fact 
that, in Cyg X-1, the gravitational power is apparently 
equally shared between the accretion flow and the jet. 
This requires discs with aspect ratio quite large, namely 
e ~ 0.1. Such solutions, like our JED ones, are nowadays 
available On the other hand, if ^ were too large the disc 
accretion rate (scaling as r^) at the inner radius during 
the hard state would be much smaller than that at the 
outer radius, while there are basically the same during 
the soft state. This would therefore lead to a ratio of 
radiative efficiencies much larger than what is currently 
observed. This upper limit on the local disc ejection 
efficiency translates into an upper limit of the jet mass 
loss, namely M jets /Ma{r out) ^ Cln(ro„t/r„,) < 50% (for 
a JED extension arbitrarily chosen to be Tout I "Tin = 100). 
At first sight, this may appear not too constraining. 
However, accretion-ejection dynamics link the jet mass 
loss to the jet asymptotic velocity. This is the reason why 
any independent measure of the jet terminal velocity is 
so important. By comparing the average velocity provided 
by the models to the values derived from observations, 
we conclude that any JED model with ^ of a few percent 
matches all observational constraints for Cygnus X-1 jets 
(velocity, power, mass loss). 

Concerning the disc itself, we found JED temperatures 
of 120-140 keV and optical depths r ~ 0.7 - 0.9 for 
Cyg X-1 parameters. These temperature/optical depth 
values are slightly larger/smaller than those generally 
deduced from ob servations, namely T^ ~ 70 — 100 keV 
and T ~ 1 — 3 (Gierlinski et al 



1997 Frontera et al. 



2001; Ibragi inovet al.[|2005[|Bel et al.||2006[|Wilms"erar 



agrees well with the observations. We thus 
expect a X-ray spectrum with the correct spectral index 
but with a cut-off (mainly controlled by the temperature) 
too high. This could be directly linked to our overly 
simple treatment of the disc radiative transfer which 
neglect for example its vertical stratification. One might 
also argue that the electron particle distribution is not 
a perfect Maxwellian. Indeed, recent computations of 
plasma radiative equilibria, including radiative and kinetic 
processes, show that a non-thermal high energ y tail in 
the particle distribu tion is generally e xpected (Belmont 
et al.[ |2008| [Malzac fc Belmont [ |2009[ ). The high energy 
photons produced by these particles could create pairs, 
then increasing the particule density and decreasing the 
temperature. These different effects will be more precisely 
studied in a fortcoming paper. 

The presence of a non-thermal high energy tail above 
the thermal Co mptonisation cut-off in the hard state 



of Cygnus X-1 (|McConnell et al 
4 (Wardzinski et al. [^ 
that the magnetic tieli 



. |2002| ) and GX 339- 
Joinet et al.l 2007) suggests 



p002[ ^^ ^ ^ 

Id is low in the Comptonizing re- 
gion. Indeed, this excess is believed to be the signature 
of a population of non-thermal electrons in the corona. 
These high energy electrons produce a self-absorbed syn- 
chrotron emission in the IR-optical bands. These soft pho- 
tons are then Comptoniscd in the hot plasma. If the 
synchrotron luminosity is too strong then it is impossi- 
ble to sustain the relatively high (~ 100 keV) tempera- 
ture of the Maxwellian electr o ns. In the case of Cygnus 
X-1, Wardzinski & Zdziarski (2001) estimated that the 
non-thermal tail implies a strongly sub-equipartition mag- 
netic field. Then based on a d etailled spectral modelling 



2006|). However, our value of the Compton parameter 



Malzac fc Belmontj (2009| and [Poutanen fc Vurm| ( |2009 ) 
infer that the ratio of magnetic to radiation energy density 
Ub/Ur < 0.3. In the framework of our model this would 
imply a magnetization that is much too small to allow for 
the production of a jet. 

However this constraint on the magnetic field holds 
only if the non-thermal tail is produced inside the ther- 
mal Comptonisation region. A possibility would be that 
the bulk of the Comptonised emission in the hard state is 
produced within the JED, while the non-thermal compton- 
isation component arises from a magnetic accretion disc 
corona above and below the standard outer disc. If so, 
the magnetic field could be strong in both JED and non- 
thermal corona. Then, during spectral state transitions, the 
non-thermal component would gradually get stronger and 
stronger as the transition radius between outer disc and 
JED moves closer to the black hole. In the soft state the 
JED disappears and the hard X-ray emission is dominated 
by a non-thermal Comptonisation in the corona. 
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Appendix A: Jet power and b parameter 

The power leaving the disc and carried away by the jets 
(from both face) is Pjets = 2 J pUpE ■ dS, where the in- 
tegration is done along one disc surface from Ri„ to Rout 
and 



E^ 



H 



*fi -a 



V 



(A.l) 



is the well known Bernoulli integral. It is the sum of all 
specific energies carried by the jet, namely kinetic, enthalpy, 
potential and magnetic. In this last term, fi* is the rotation 
speed of the magnetic surface and rj is the mass flux per 
magnetic flux unit. 

In the ideal MHD regime realized in jets, E is an in- 
variant that can be evaluated at the disc surface. In a thin 

2 

disc (e <C 1), it writes E ~ ^(2Abp — 3), where Xgp is 



the Blandford & Payne ( 1982 ) magnetic lever arm param- 
eter. For the models considered here, Xgp varies between 
10 a nd 10 so that the magnetic contribution dominates in 
Eq. (A.l). The jet power leaving the disc is therefore mostly 
in the form of an MHD Poynting flux. 



P. 



Pa 



jets — ^MHD 



2 / n,R^ 2TTRdR 



(A.2) 



which is exactly Eq. (13 1. The toroidal magnetic field at 
the disc surface scales as B'f = —FiHoJoH where Jo is the 
radial electric current density at the disc midplane, with 
Fi ~ 0.4—0.6 for most published solutions. In a JED, where 
the dominant torque is due to the jets, the disc angular 
momentum conservation is 



J „ PoUo 5 2 1 o 

Jo^z — P Jui ° ^ -^PoUo^l-o 



(A.3) 



where the subscript "o" stands for values evaluated at 

the disc midplane. Writing Ma{R) = —AnR / pu,.dz — 

Jo 
4ttRpoUoHF2, with F2 a number of order unity found to 

vary between 0.4-0.6 also for most of our solutions, we ob- 
tain 



jets 



R^n 



"^ ^]^ J jet 



(A.4) 



In this expression, fjet{R) = — 7;; — 7;; — is a priori a func- 

tion of the radius. However it cannot vary much if ejection 
is to occur on a large radial extent. Under this ass ump tion, 
our self-similar scalings apply and we obtain Eq. ( 14 ) with 



F1/F2 ii* ilo ''* "O 



(A.5) 



While this quantity is nearly independent of ^, it is clearly 
a function of the disc thickness, hence its temperature. The 
rotation of the disc material at the disc midplane scales as 

ilo = ^K\/^ — fe^ — PfJ-^, where the second term is due to 
the radial pressure gradient and the third to the magnetic 
tension (p ^ B^ /B^ ^ 1). Thus, the thicker the disc, the 
larger the deviation from Keplerian and the smaller b. 

But the dominant effect is due to the rotation of the 
magnetic surfaces fJ* . Indeed, fi* is defined only in the up- 
per layers of the disc once ideal MHD becomes valid. The 



disc density has therefore decreased significantly, enhancing 
thereby the outward effect of the magnetic radial tension. 
This tends to decrease the angular velocity fi+ of the mat- 
ter and thus of the magnetic s urfaces fi^, (see for instanc e 
the vertical profiles shown in Ferreira & Pelletier 1995). 
This behaviour can be also understood with the magnetic 
torque. This torque is first negative before becoming posi- 
tive at the disc upper layers. The disc is therefore first spun 
down and fi decreases. 

Numerical solutions obtained with isothermal magnetic 
surfaces provide h ranging from 0.99 for thin discs with 
£ smaller than ^ 10^^ to roughly 0.5 for slim discs with 

£-0.1. 
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Fig. A.l. Parameter space for the ejection efficiency ^, the 
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Mach number mg for the accretion-ejection solutions of 



Ferreira (19971 used in this paper. The fines correspond 
to different values of the disc aspect ratio H/R: solid lines: 
H/R = 0.1, dashed lines: H/R = 0.01, dot-dashed lines: 
H/R = 0.001. For this set of solutions, the turbulence level 
parameter a„i is taken equal to unity. 
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